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The role of fluctuation-dissipation relations (theorems) 
for the magnetization dynamics with Landau-Lifshitz-Gilbert 
and Bloch-Bloembergen damping terms are discussed. We 
demonstrate that the use of the Callen-Welton fluctuation- 
dissipation theorem that was proven for Hamiltonian systems 
can give an inconsistent result for magnetic systems with dis- 
sipation. 
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I. INTRODUCTION 

The study of linear stochastic magnetization dynamics 
is of great importance in applications to nano-magnetic 
devices and ultra-thin films. In general, gyromagnetic 
magnetization motion around an effective field is ran- 
domly forced by fluctuations on spins by means of inter- 
action with a thermal bath (phonons, magnons, conduc- 
tion electrons, impurities, etc.). Fluctuation-dissipation 
relations give a very useful tool to find a correspondence 
between the dynamic variable fluctuations, temperature 
and dissipation for a given magnetic dynamic system 
(see, e.g., [1], [2], [3], [4], [5]). 

The most frequently used fluctuation-dissipation re- 
lations are known as the Callen-Welton fluctuation- 
dissipation theorem (FDT) [6], [7]. There are two in- 
tegrals that express this theorem. The first one gives the 
correlation function of the dynamic variables in terms of 
an integral of the system susceptibility. The second form 
relates the correlation function of the applied noise to 
an integral of the inverse susceptibility. This theorem 
is proved for rather general assumptions (we will discuss 
them later) and therefore appears to be very attractive 
for general problems. The aim of this paper is to discuss 
the inapplicability of the second relation to magnetiza- 
tion dynamics in ferromagnetic systems. In particular, 
we demonstrate that Callen-Welton FDT gives mathe- 
matical inconsistency and can not provide an argument 
to discriminate between different damping models. 



II. FLUCTUATION-DISSIPATION THEOREM 

The Callen-Welton fluctuation-dissipation theorem is 
proved under very general assumptions for Hamiltonian 
systems of an arbitrary type (see, e.g., [8], [9]). These 



systems have no internal dissipation at all. The energy 
V of the external perturbation acting on the system is 
taken in the form: 



v = -5>i(*)/j(*). 



(i) 



where Xj (t) is the j-th component of the dynamic variable 
and fj (t) is the j-th component of the external alternat- 
ing field. All dynamic variables are assumed to be equal 
to zero in the absence of the perturbation and the linear 
responses are defined by 



z 

x A t )=Y. I K jk (t-T)f k (r)dr, 



(2) 



where Kjk(t) is a memory function that depends on the 
properties of the dynamic system [10]. A "thermal bath" 
in this theorem is modeled as a set of periodic external 
fields, which are absorbed by the dynamic system. On 
the other hand, the external fields stimulate radiation 
from the system and cause a loss of energy. For a periodic 
external field 



/i(*) = 5 [/oje- <w *+/5 J .e & -*] 



one can rewrite Eq. (2) in the form: 



(3) 



(4) 



where 



Xjk{u) = J Kjk(t) exp(iujt)dt 




(5) 



is the susceptibility tensor. 

The change of the total energy £ (which includes the 
perturbation energy V) in the dynamic system that plays 
an important part in this theory is expressed as 



d£ 
~dl 
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dfjjt) 
dt 



(6) 



Averaging this equation over the period of the external 
fields (3) and taking into account Eq. (4), one obtains 
the following expression for the dissipated power: 
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Q = TE[^M-Xfc J -M]/o J -/oV (7) 

In the derivation both the absorbed and radiated aver- 
aged powers are usually expressed quantum mechanically 
in terms of transitions probabilities (per unit time) be- 
tween infinitely narrow energy levels, i.e., no linewidth 
or damping is taken into account. Temperature is intro- 
duced by the thermodynamic weights of the energy lev- 
els. This gives a thermal averaging over a thermal bath. 
Summing up over frequencies ui, one can obtain the bal- 
ance between absorption and irradiation. This balance 
expressed in terms of dynamic susceptibility, gives the 
Callen-Welton FDT. In the classical limit it can be writ- 
ten as: 



{ Xj (r)x k (0)) 



k B T 
2tt 



f 
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Here ks is the Boltzman constant, and T is the tempera- 
ture. (. . .) denotes both thermal averaging and averaging 
over random phases in the periodic fields. 

A second form of the Callcn-Wclton FDT is also intro- 
duced: 



(fj(r)fk(0)} 



k B T ' x 

~2^r 



(9) 



This form results from the reversibility of the linear re- 
lation 



to 



(10) 



(11) 



It should be noted that the first (8) and second (9) rela- 
tions can be interpreted differently Namely in the first 
case we have a reasonable correspondence between the 
correlations of dynamic variables and the susceptibility 
of the dynamic system. Changing the dynamic system 
properties, and hence the susceptibility we change the 
correlations of the dynamic variable. A specific external 
noise mechanism is excluded and the role of the thermal 
bath is included implicitly by assuming a system in ther- 
mal equilibrium. This gives a linear dependence on the 
temperature T. 

On the other hand, in the second case the noise correla- 
tions are associated with the inverse susceptibility of the 
dynamic system. In other words, two principally different 
physical characteristics, the noise, which results from a 
thermal bath dynamics and the susceptibility, which de- 
scribes just the properties of the dynamic system, should 
correspond to each other. From Eq. (9) it follows that 
by changing properties of the dynamic system, one can 
change the noise correlations. 



However, physically it seems impossible in general that 
the dynamic susceptibility determines the noise variance. 
It is not clear that the correspondence between (10) and 
(11) applies for general random processes. This follows 
from the fact that the integral (2) depends on the form 
of stochastic integration because /fc(r) are random vari- 
ables. For example, the Ito and Stratonovich approaches 
to stochastic integration have no general relationship be- 
tween each other in the case of multiplicative noise [11]. 

In applications to systems with dissipation, the Callen- 
Welton fluctuation-dissipation theorem may be used just 
as an approximation, which must be validated. For some 
dissipative systems the Callen-Welton FDT gives rea- 
sonable estimates. It covers, for example, the Einstein 
relation for Brownian motion and the classical Nyquist 
formula for voltage fluctuations. However, these particu- 
lar cases do not prove that the theorem is applicable for 
any linear system with dissipation. Some criticism of the 
Callen-Welton FDT is given by Klimontovich [9], [12]. 



III. FDT AND MAGNETIC DAMPING 

A commonly used theoretical approach to magnetic dy- 
namics, which is purely phenomenological, is based on 
the Landau-Lifshitz equation [13], or its modification in 
the Gilbert form [14]. Local and isotropic phenomenolog- 
ical damping terms and corresponding local random fields 
are assumed to describe a "thermal bath" in a magnetic 
material. On the other hand, a microscopic approach pre- 
dicts non-local field fluctuations [15] and magnetization 
dynamics in the form of the Bloch-Bloembergen equa- 
tions: the Landau-Lifshitz-Gilbert equations occur only 
for uniaxial symmetry [16], [17]. 

Recently Smith [18] using the second form of the 
Callen-Welton theorem (9) claimed that fluctuation- 
dissipation relations provide a means for discriminating 
between alternative phenomenological magnetic damping 
models. Here we demonstrate that, by using Eq. (9), the 
fluctuation-dissipation relations for stochastic dynam- 
ics with Landau-Lifshitz-Gilbert damping and Bloch- 
Bloembergen damping are inconsistent. 

For simplicity we shall consider a single-domain parti- 
cle with the magnetic energy £ in the vicinity of equilib- 
rium: 



£ 



M S V 



~ir m l + -ir m l ~ h x m x - h v my. (12) 



Here m x and m v are the transverse components of the 
unit vector m oriented along the magnetization, M s is 
the saturation magnetization, V is the volume of the sam- 
ple, H x and H y are the Kittel stiffness fields (in general 
H x 7^ H y , for example, in a thin film). For a general 
stochastic problem h x (t) and h y (t) are random variables 
arising from the interaction with a thermal bath. To 
illustrate the utilization of Callen-Welton FDT one as- 
sumes that h x (t) and h y (t) are equivalent to the external 
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alternating fields fj in Eq.(l) that leads to the forms of 
Eqs.(8) and (9). 

The linearized Landau-Lifshitz-Gilbert equation can 
be written in the form: 

i (m x +a m^j = -H y m y + h y (t), 
■^(m y -am x ^J = H x m x - h x (t), (13) 
or, equivalently as 

lm x = -aH x m x -H y m y + h y{ t) + ah x (t), 

- m y = -aH y m y + H x m x - h x (t) + ah y (t), (14) 



where 7 = 7/(1 + a 2 ), 7 is the gyromagnetic ratio and 
a is a dimensionless damping parameter. From Eq. (13) 
one can obtain the inverse susceptibility 
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lM s V 



-iau> + jH x 



-iau) + -fH y 



(15) 



where hj{oj) = xjk (u)mk(uj)M s V; mk{w)M s V is the k- 
th component of the magnetic moment. Substituting Eq. 
(15) and its Hermitian conjugate into Eq. (9), the second 
form of Callcn-Welton FDT can be expressed as 



(h x (T)h x (0)) (h x (T)h y (0)) \ _ 2k B Ta ( 1 \ 
(h y (r)h x (0)) (h y (r)h v (0)) J ~ jMs V [o 1 ) ^ 

(16) 



The case of linearized Bloch-Bloembergen magnetiza- 
tion dynamics is described by 



- m x = \—m x - Hymy + h y (t), 

7 7 J 2 



1 

7 



1 



- m y = —m y + H x m x - h x (t), 

I-L2 



(17) 



where T 2 is the relaxation time. From this equation we 
can obtain the following inverse susceptibility: 
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lH x 



iuj — T 9 1 



X]k{Uj) JM S V {-iuj + Tz 1 1 H y 
Application of Eq. (9) yields FDT: 



(18) 



the Bloch-Blocmbcrgcn form of damping is not physical 
due to the non-diagonal form of (19). However, as we 
show below, there is an inconsistency in the use of the 
second form of the Callcn-Welton FDT. 

Let us now consider the uniaxial case when both dy- 
namics should coincide. We examine the dynamics of a 
soft micromagnetic particle (no anisotropy) in an exter- 
nal magnetic field Ho and small damping approximation 
q«1 (the most interesting case) , when terms ~ a 2 can 
be neglected. Thus H x = H y = H , and the Landau- 
Lifshitz-Gilbert equation (14) is reduced to: 



(Mt)MO)) (hx(r)hy(0)) \ k B T ( -1 \ sgn(r) 
(h y (r)h x (0)) (h y (T)h y (0)) J lMs V \ 1 J T 2 

(19) 



Note that the difference between Eq. (16) and (19) lies 
in the tensor on the right hand side of each equation. 
Because of the diagonal tensor, Eq. (16) implies no cor- 
relations between h x and h y . On the other hand, Eq. 
(19) does show such correlations. One can argue that 



1 dm x 
7 dt 
1 dm 



-aH m x - H m y + h y (t) + ah x (t), 



^ = -aH m y + H m x - h x (t) +ah y (t). (20) 



7 dt 

So far as h x (t) and h y (t) represent two independent ran- 
dom fields, their linear combinations are also random. 
We can consider the following orthogonal transformation: 



hy(t) = 



h y (t) + ah x (t) 



h y (t) + ah x (t), 



_ Mt) = + ^ - hx[t) + ah y {t). (21) 

V 1 + ot A 

Thus, random fields h x (t) and h y (t) are independent and 
Eq. (20) can be rewritten in the form: 



1 dm x 
7 dt 



= aH m x - H Q m y + h y (t), 



1 dm 



v _ 



7 dt 



= aH m y + H m x - h x (t). 



(22) 



From Eqs. (16) and (21) one can easily calculate pair 
correlations for h x (t) and h y (t) (the terms <~ a 2 are ne- 
glected): 



(MT)MO)) (h X (T)hy(0)) 

(h y (T)h x (0)) (h y (r)h y (0)) 



2kTa 



1 
1 



S(r). 
(23) 



Note that the Eq. (22) is an exact mathematical analog 
of the Bloch-Bloembergen equation (17) (l/T 2 = ajH ) 
for the transverse magnetic components with random 
fields h x (t) and h y (t). Thus, we see that the fluctuation- 
dissipation relations (23) and (19) are principally differ- 
ent in the case when they must coincide (Gilbert and 
Bloch-Bloembergen dynamics are the same) . The second 
form of the Callen-Welton FDT can not be used to verify 
the validity of one form of damping versus another. 

The origin of the inconsistency is in the use of the sec- 
ond form (9) of Callen-Welton FDT, applied to systems 
with dissipation (Landau-Lifshitz-Gilbert and Bloch- 
Bloembergen equations). Without dissipation (a = 0) 
this inconsistency disappears. It should be pointed out 
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that Eqs. (20) and (22) describe stochastic magneti- 
zation dynamics with (h x (t)) = (h y (t)) = (m x (t)) = 
(m y (t)) = in accordance with thermodynamics. 

Note that for no damping the dynamic susceptibility of 
the Landau-Lifshitz-Gilbert equation (15) and the Bloch- 
Bloembergen equation (18) are identical in form. In the 
stochastic case with random fields h x (t) and (h y (t)) at 
first glance Eqs. (14) and (17) appear to be different. 
However, we have shown that by a simple transformation 
(21) both equations are equivalent in the case of uniaxial 
anisotropy. 

Application of the first FDT (8) to both Landau- 
Lifshitz-Gilbert and Bloch-Blocmbcrgen stochastic dy- 
namics does not give such a strong inconsistency as does 
the second form. Using (8) with (15), we obtain the fol- 
lowing averages for the Gilbert dynamics: 



(to x (t)to x (0)) (m x (T)m y (0)} 
(m y (T)m x (0)) K(T)m„(0)) 



(24) 



nM s V 



^gHI 2 

w 2 (l + a 2 ) + ( 7 H y ) 2 -iun(H x 
ioj-f(H x +H y ) u 2 (l + a 2 )- 



^H y ) 
{lH x f 



where Dq{u) = loq — uj (1 + or) — iauj"f(H x + H y ) and 
uJq = -f 2 H x H y is the ferromagnetic resonance frequency. 
On the other hand, the use of (8) with (18) gives: 



(m x (r)m x (0)) (m x (T)m y (0)) 
(m y (T)m x (0)) (m y (T)m y (0)} 



(25) 



lT^k B T 



ttM s V 



f 



du> 



i2uj"fH y 



iuj\D BB {uj)\ 2 

uj 2 + uj 2 + T- 2 
r^ 2 i2w-)H x 



where D bb (oj) = uj 2 — ui 2 + T^ 2 - UuT^ 1 . 

In general, magnetic correlations (24) and (25) differ 
from each other. However, in the most important case of 
the noise power (r = 0) both (24) and (25) are reduced 
to 



K(0)> 

(771^(0)771,(0)) 



(777,(0)777^(0)) 

K(o)> 



k B T ( H- 1 
M H V V 





(26) 



This equation is completely consistent with thermody- 
namics, namely, with energy equipartion: 



M„VH X 



{m 2 x ) 



(27) 



been derived by Kubo [19] for the permeability of a dy- 
namic system and by White [2] for the susceptibility of 
a general magnetic system. A standard method to study 
stochastic dynamics in systems with dissipation is the 
Langevin approach [9] , [20] . Application of the Langevin 
approach utilizing specific dissipation mechanisms can be 
found in [15], [16], [17]. 

The analysis of damping in magnetic systems has had 
a long history (e.g., [21], [22]). Each spin wave (magnon) 
interacts with a so-called, thermal bath, which consists of 
magnons, phonons, conduction electrons, impurities, etc. 
Analyzing the processes of relaxation, one can find the 
spin-wave damping (see, e.g., [21]) and the correspond- 
ing thermal noise. The microscopic stochastic differential 
equation (SDE) is shown to be of the Langevin form for 
a damped harmonic oscillator for a wide class of relax- 
ation processes [17]. Note that the collective description 
works even in the case of local interactions. A local defect 
or impurity perturbs the band structure of the magnetic 
crystal (see, e.g., [2]) and affects the spin- wave spectrum, 
whose imaginary part is responsible for the damping of 
collective magnetic excitations. 

The occurrence of delocalized damping has been di- 
rectly demonstrated in the problem of two coupled mag- 
netic grains with local thermal baths [15]. Stochastic 
forces are uncorrelated in the spin- wave coordinates, but 
become correlated when the SDE's are expressed in the 
original magnetization coordinates. Because of the sys- 
tem interactions, even though the coupling to the ther- 
mal bath may be purely local, there is no physical re- 
quirement that the stochastic forces in the magnetization 
coordinates be uncorrelated. 

In summary, we have argued that the second form of 
the Callen-Welton fluctuation dissipation theorem does 
not correctly apply to damped systems, has inconsisten- 
cies in its application and thus can not distinguish be- 
tween relaxation models. 
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IV. DISCUSSION 

Other forms of the fluctuation-dissipation relations, 
which arc similar to the first Callen-Welton FDT, have 
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